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D I Abstract 

I The purpose of this note is to point out ambiguities that appear in the calculation of angular 

, momentum and its radiated counterpart when some simple formulae are used to compute them. 

We illustrate, in two simple different examples, how incorrect results can be obtained with them. 
Additionally, we discuss the magnitude of possible errors in well known situations. 

I ! 1 Introduction 
^ . 

It is well known that the notion of total angular momentum in general relativity is sensitive to the so 
called problem of supertranslation ambiguities. This topic involves subtleties and complications inherent 
in the infinite dimensional nature of the asymptotic symmetry group of asymptotically flat spacetimes, 
namely, the BMS group. On the other hand, there is a clear need to compute these quantities to extract 
' valuable physical information in numerical constructed spacetimes to connect with observations of a 

\^ ■ variety of sources. A few proposals exist to do this either locally [U [21 El HI E] or asymptotically [H El [71 

O; [Hl[a[ini[IIl[I20. The dependence of any notion of relativistic angular momentum on supertranslations 

has been the matter of numerous studies. The problem is complex and, unfortunately, most definitions 
of angular momentum either show supertranslation ambiguities or require of additional structure (at the 
OO . expense of reducing generality in their application) . 

In spite of all these difficulties, the need for tackling concrete calculations of angular momentum 
in numerically constructed spacetimes has motivated the use of simple and practical formulae for the 
computation of angular momentum. These formulae lend themselves for a straightforward computation 
and, in relatively simple scenarios, have given the expected answers. One option that has been employed 
' in several works has been suggested in |13| , where it is proposed to compute the angular momentum from 

the expression 
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= Y^SRe dS<i>[,] [2*0 - 2a°da° - g(a°a°)] | ; (1) 

which is a truncated version of the one defined in [8J; in fact, a supermomentum part is missing. As 
mentioned above, supertranslations issues can obscure angular momentum definitions (and expressions 
derived from them) and indeed, as we show in the next sections, eq. ((T]) is susceptible to them as well. 
We discuss in the next examples two simple cases where expression ([1]) gives an incorrect result. 
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2 Examples 



2.1 A supertranslated and boosted Schwarzschild 

Let us consider, for simplicity, the case of a spacetime describing two orbiting black holes which after 
merger give rise to a non-spinning black hole. Furthermore, we consider a case where the final black 
hole has a kick. This scenario, can easily be realised with a non-equal mass binary with individual spins 
adjusted such that the final spin is zerqj. Therefore, the black hole will generically appear in the grid 
which is supertranslated from that where the black hole would be at rest. 

To simplify the analysis, we further assume that a natural coordinate system can be constructed such 
that, asymptotically, the system can be surrounded by round spheres; that is, one can construct a Bondi 
type coordinate system from it. The delicate issue is that this inertial asymptotic coordinate system 
results, generically, supertranslated and boosted with respect to the intrinsic asymptotic coordinate 
system which has the spherical symmetry of the remaining black hole. This in turn, implies that the 
asymptotic fields acquire a non-trivial dependence on the supertranslation and boost. To compare what 
would be obtained in both systems, we denote with primes the asymptotic coordinate system with 
spherical symmetry and unprimed that naturally constructed in the numerical implementation; then one 
would have the BMS transformation 

u = KiC, C)iu' - 7(C', C')) + ^^i^^^^l^ (2) 



where K is the Lorentz factor associated with the boost, 7 determines the supertranslation and (a, 5, c, d) 
are complex constants which determine the boost; for details see |19) . 
In what follows it will be convenient to use the notation 

Although the calculation of expression H]) on a spherically symmetric section gives the expected null 
result we will show next that it gives a non-zero result when calculated at a supertranslated and boosted 
section. 

As explained elsewhere |20[ 121] the above transformation implies a transformation of the null frames, 
and in particular, the components of the Weyl tensor and the asymptotic shear are affected. More 
concretely, one has 




if3 \ K ^ \ K 

Only ^'2'^ = —m is different from zero in the spherically symmetric frame; therefore one has 

*? = ge'^-g'^; (7) 
where m is the mass of the resulting black hole. Similarly the shear is 

^--e'^'d'^j. (8) 
K 

Let us now calculate expression ([I} in detail, at a u =const. section, and discuss what is obtained. 
First, note that the factors in the integrand, are spin weighted quantities with spin weight -1. 
Therefore they can be expressed in terms of an edth bar operator acting on a spin weight quantity. 



straightforward calculation can be employed to adjust the physical parameters so that this scenario is realized as 
discussed by the simple model |14j or the fitted expansions presented in [151116) . The resulting kick velocity can be estimated 
by suitable fits as in [T5l[T71[T8l . 
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namely = Sa^jj ; where, to pick the angular momentum components ajjj should be imaginarv|ll|. 
Consequently, each a^jj is a quantity of spin weight satisfying SSa^jj — —ay^y 

We should now express the integration appearing in ([T]) in terms of the primed, spherically symmetric 
Bondi system. Let us observe that for all these quantities a[i] , due to the transformation properties of 
the edth operator [22], one has 



8'a'; (9) 



where a' ~ a'((^' , (.') is understood in terms of the primed angular coordinates. Next, one must also take 
into account that dS = K^dS', then the first term of ([T]), involving "^fl, is 

5Re = 5Re ^$|,] ^e^' 3' [K{u' - j)]K^dS'^ . (10) 

In order to simplify the discussion we can assume that the calculation is performed at the section u = 0, 
which is equivalent to u' = 7; therefore, one has 9' [K{u' — 7)] = {u' — 7)9' K — Kd' 7 = —Kd' 7; which 
implies 

(11) 



Since the factor K~^d'j has spin weight 1, one again can think of a spin-weight potential A such that 
Q' A = K~^B'^. Therefore, if the quantity A were real, this term would not have any contribution, since 
one would have to evaluate 



3fie |£ ^[.j^'JdS'l = - 5Re jsm £ B'a[^&AdS'^ = 5Re jsm £ d'B'a[.^^AdS' 



^elSmj) a[^]AdS' > ^ 0; 



(12) 



since aj-j's are imaginary and we are assuming at the moment that A is real. Unfortunately, it is not 
difhcult to prove that actually A has an imaginary part; and furthermore, this imaginary part is only 
zero when the gradients of -j^ and 7 are proportional; something that is not true in our case. 

To see this, we next study the equation BA = v'^B^; where v = K~^, and assuming A is a real quantity 
we will reach a contradiction. Notice that for a real spin weight zero field A, Im [BSA] = 0; since §3 is a 
real operator when acting on a spin weight zero quantity. This last relation can be written in terms of v 
and 7 as 

Bv^B-f - d-/dv^ = 0. (13) 

Let us denote the gradients on the sphere of the quantities and 7 by VaW^ and V67. Then the 'cross 
product' of these vectorial quantities is defined by e"''VoW^Vb7; where e"'' is the surface element of the 
unit sphere. At this point it is convenient to express the gradients in terms of the edth operator; for 
example for v"^ one has 

Vav^ = ~{dv^)ma - idv^)ma; (14) 
and similarly for 7. Then the cross product is 

e''^VaW^Vf,7 = {Sv^d-f ~ d-ydv^) t^^mafru,; (15) 

which should be zero if A were real. But this is true only if the gradients are proportional. Consequently 
the term ()lip has a non-zero contribution to the quantity J[i] . 

We now turn our attention to the second term in expression ([1]) which involves quadratic terms in 
a. One could also in this case think of a potential X such that cr" = Then, one can prove that 

if X were real, then the second term of ((!]) involving the cr's would not contribute to the integral[2^. 
However, following a similar line of arguments as above, one can prove that X will have an imaginary 
part if 

--ie^^Vl, (16) 

-recall this is the factor appearing in ([S])-. 

Therefore, this second term will also have a non-zero contribution to the quantity J^^j; and since it 
includes independent information not contained in the first term it can not cancel the first contribution. 
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Let us discuss further what we have found up to this point. If the adopted coordinate system, for 
instance the one most naturaUy constructed from a numerical simulation, results supertranslated and 
boosted with respect to the final frame adapted to the spherically symmetric final black hole an incorrect 
result will be obtained for the angular momentum of the system. Notice that this conclusion is more 
subtle than the obvious 'Newtonian'-type observation stating that the angular momentum calculated at 
frames translated with respect to each other are different. In fact, even if the coordinate system is only 
boosted and supertranslated -without a translation-, with respect to the one adapted to the spherically 
symmetric final black hole then J[i] turns out to be different from zero in the Schwarzschild geometry. 
This is to be contrasted with what would be the case within special relativity paradigm where the absence 
of a translation does not affect the angular momentum calculation since the Lorentz transformation of a 
zero tensor is a zero tensor. 

2.2 A radiating second order solution 

As a second example let us consider the spacetime given by the metric 

(V m\ 
-2—r + Kv-2—\du^ + 2 du dr - —r-^dC, dC; (17) 
V r I ^ Pq 

where a dot denotes a derivative with respect to u, Ky is the Gaussian curvature, given by 

Kv ^2VWV -2dV BV + V^, (18) 



of the 2-metric 

^ dC dC; (19) 

and Po{C, C) is the conformal factor of the unit sphere. This metric was used by Robinson and Trautman 
to characterize spacetimes satisfying Einstein equations admitting a congruence of null geodesic without 
shear and twist [21]. 

The scalar V{u, Q can be expressed in the form V = 1 + A(u, (, (); so that A = would correspond 
to the Schwarzschild spacetime. The function A is required to satisfy the equation 

-3 M A = {l + 4A)dWA~B^Ad^A. (20) 

This spacetime is known as a Robinson- Trautman geometry [25]. By requiring (|20p one obtains a 
solution of vacuum Einstein equations in second order in A. The Robinson- Trautman spacetimes are 
understood as representing a central object which has been perturbed and are known to decay in time to 
the Schwarzschild geometry. 

Suppose one has obtained this spacetime and would like to calculate the quantity expressed in ([T]) at 
the section u — uq- Then, since the asymptotic coordinate system (u, C, C) is not inertial (not Bondi); one 
should look for the transformation to a Bondi frame. Let us choose a Bondi coordinate system {u, ^, 
such that u = coincide with the section u = uq. With this choice one has that on this particular section 
the Bondi shear is zero; = 0[1T]. Furthermore, since on this section the null frames are proportional, 
one also has VP^ = 0. Therefore in this section one obtains Jj^j = 0. 

However in any other Bondi section given by m 7^ one will obtain J[j] ^ 0. This can be seen due 
to two reasons. First one can notice that the shear with respect to the Bondi system is not stationary; 
in fact one can prove that tT^ = (S^y)F~^[25]. One can convince oneself that a can not be expressed 
in terms of a real potential X , that is, a = d^X as studied in the previous section. In other words such 
an X has an imaginary term. Therefore one deduces that there is always a flux, which in general will 
include the time variation of J[i], depending on the arbitrary initial data V. On the other hand one can 
also infer that ^'j' will also be different from zero, at any other Bondi section; due again to the existence 
of radiation. 

Notice that J[i] 7^ at any other Bondi section u ^ can also be deduced by taking (Bondi) time 
derivatives of J[j]. One can see that at m = 0, one has Jj^j - = 0, Jj^j - - = but J[j] ^ 0. 

The weaknesses of formula (1) to calculate the angular momentum also give rise to odd results. Notice 
that the asymptotic value of Jj^j , for -u 00 is zero; since in this regime, u will be supertranslated with 
respect to the (w, C, C) frame; but not boosted (see the previous section). So, in this case one would have 
the curious situation that the initial value of J[i] is zero, a non-zero flux to the future of the initial section 
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and an ending with a zero value for J[j] in the asymptotic regime. This alone makes it difhcult to identify 
J[,j] with a physical notion of intrinsic angular momentum. 

One instead could use the Bondi system (u*, C*j C*) which is defined so as to coincide asymptotically 
with the R-T system {u,(^,() in the regime u — s- oo. Then one could calculate the quantity Jj^j for any 
u* ^constant. In this case one would obtain that in the limit u* — s- oo, Jj^j would vanish. However for 
any finite u* =const., the evaluation of Jj^j would give a non zero value. 

This again stresses the difficulties that one encounters in trying to use the expression ([T]) to obtain 
intrinsic information of the central object. 

3 Estimate of the error in spin calculations 

In order to provide a diff'erent perspective to the quantity J[j] , in the next subsection we recall an 
unambiguous definition of total intrinsic angular momentum at future null infinity. Actually since this 
is the only definition we know of that solves simultaneously the supertranslations ambiguities for center 
of mass and intrinsic angular momentum [11 j -without assuming further structure-; we will use it as 
reference for comparison. 

3.1 An unambiguous definition of intrinsic angular momentum 

Arguably one of the most valuable physical information aimed to be extracted from a given spacetime is 
that of intrinsic angular momentum as it relates directly with observable signatures at both gravitational 
and clcctromagnetical signals produced in a number of systems. Tied to defining such quantity is the 
problem of defining the center of mass of the system with respect to which the intrinsic angular momentum 
can be defined; consequently, the definition of both concepts must come together. A solution of these 
difficulties has been formulated in [TT]; where the expression to calculate the charges associated to BMS 
generators is given by, 

Os™H =^el^8j^ (- W2 + 2aoddo + 9 {a„<Jo)) + 2wi (*° + ao(fo + S^-q)) dS'^ ; (21) 

where wi and W2 are components of the two form w that is determined by the particular generator of 
BMS transformations. As with the expression appearing in [H], one can use this type of expression to 
calculate the total momentum, total supermomentum and total angular momentum. However, the notion 
of intrinsic angular momentum can only be obtained if one calculates the angular momentum on the 
center of mass sections, here denoted by Scm- 

If one calculates the angular momentum in a radiating spacetime at any other section, one would 
obtain a quantity with angular momentum reminiscence but with unclear intrinsic physical meaning. 

3.2 How important are in practice these considerations in the computation 
of angular momentum? 

We want to estimate now the errors in the value of angular momentum by using an unfortunate choice of 
frame and section. To do so, we want to have an estimate of how supertranslated are the "center of mass" 
sections with respect to the sections that are adapted to the coordinates used in numerical computations. 

The astrophysical model that we have in mind is similar to the system that was presented in subsection 
12. H i.e, we have some compact objects (for example a binary system) that are undergoing a merge, and 
where the coordinate grid is in some way following the system until the moment of merger. After such 
event, in general the system will end up with a unique boosted (quasi)-stationary compact object (in 
the case of black hole collision, we will finish with something similar to a boosted Schwarzschild or Kerr 
black hole). Then, the new "center of mass" sections will be supertranslated and boosted with respect 
to the originals (those adapted to the numerical grid). 

Sections that characterise rest frames are known as 'nice' sections |26). The prescription to find rest 
frame sections uses a particular notion of supermomentum. 

The supermomentum at S is defined pS] by 

Pi,r,{S) = YirniCX) ^ dS^; (22) 

s 
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where ^ = 'I'2 + oa + 9^cr and Yim are spherical harmonics. 

Given another section S of future null infinity, one can find another Bondi system (u, C) such that 
S is determined by u = 0. The relation between the new and the original Bondi system is given by a 
BMS transformation, as given by 0, © and (|3]); where now we use tildes for the new coordinates and 
un-tilde for the original coordinates. Note that S can also be determined by u = 7(C, Q . 

The section 5* is said to be of type nice 26J if all the 'spacelike' components of the supermomentum, 
when calculated with respect to the adapted Bondi system, are zero; 

nm(5) = for /^O; (23) 

and therefore the only non- vanishing one, Pqq(S)^ coincides with the total Bondi mass at S. 
One can prove that ^ transforms under a BMS transformation as 

* = (VI/ - g^e^) . (24) 

Then, for a section to correspond to a 'nice' section it must obey ^ — constant as otherwise some moments 
of Pjm, for I 7^ 0, will be non-zero. The nice section equation can be understood as a condition for 7(^, Q 
and K. 

It was indicated in ref. [26j that equation ([24]) can be expressed by 

g2g2^ ^ r c:)dv; + ^ 0, C, C) + ^(7; C, lfM{^). (25) 

Jo 

where Af(7) is the total mass at the section u = 7 and if (7; Ci C) is the conformal factor of the BMS 
transformation that aligns its timelike generator (defined as poo — y"nn^i26J) with the Bondi momentum 
atu = 7(C,C)-. 

Note that — (T^ct", and then it is proportional to the content of gravitational radiation. Now, let 
us introduce a parameter A to measure this content of radiation, i.e., let the integrated flux, appearing 
in the first term of the right hand side of ([25]) be of order A. 

Suppose then, that we start by describing our compact bodies system with sections that are close 
to 'nice' sections. Then, after the merger, a new compact object is obtained, and due to the emission 
of gravitational radiation, the new 'nice' sections that 'follow' the system will, in general, not only be 
boosted but also supertranslated by a quantity 7. Now, let v be the kick velocity of the final compact 
object; them, one can observe that k, 1 + O (u) . On the other hand, by computing the flux of the 
Bondi momentum expression we can see that 

-o(A); (26, 

and *(m = 0, C, C) ~ "^^^ then 

92927 w 0(A) + 0(yM) w 0(yM). (27) 

From these estimates, and recalling that w 0(M), and using eq.([7]), we have that A^*" w ©(uM^). 

This will be the order of the error that one would have in the computation of angular momentum 
by not adopting the correct sections. We conclude then, that the errors will be, at least, as large as the 
relativistic velocities of the astrophysical system, i.e., AJ[j] k, 0{vM^). Consequently one can understand 
why the use of the flux of eq. ([T]) , as calculated in [13] , employed in several works [571 UHl I2S 1301 131] 
(where the velocities of the astrophysical systems under consideration are low, and therefore the error 
estimate AJj^j is small) can produce sensible results. However one should be cautious with the use of 
formulas like ([T]) in cases where the individual objects move relativistically (as in the studies of high speed 
collisions [32]), large production of gravitational radiation or producing an object with a small final spin 
(where the errors described above will be important). 

Finally, we stress that all these considerations are true if one assumes that the numerical grid and 
gauge gives rise to an extraction frame/coordinate system which is of Bondi type (which in general will 
not be the case). If it is not of this type, one must consider extra gauge effects like those reported in (20j 
and |21j in the computation of radiation flux and total momentum respectively. 
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4 Final comments 



The difficulties shown above are related to the nature of the notion of angular momentum. In special 
relativity the angular momentum is expressed in terms of an antisymmetric tensor 

jab ^ gab ^ paj^b _ pa pb . (28) 

where P° is the total momentum, 5"°^ is orthogonal to the momentum and i?° shows the dependence of 
the angular momentum with respect to some reference point, or origin. 

The tensor S°'^ does not depend on the choice of origin and therefore is the intrinsic angular momen- 
tum. At the center of mass line one has J"'' = S'"'' and the original expression of the angular momentum 
gives the intrinsic value. 

The expression ([T]) does not have this translation/supertranslation dependence; and therefore its 
usage is ambiguous since there is no prescription of what a center of mass is. Instead, we have shown 
that one could use a definition of angular momentum that does show this translation/supertranslation 
dependence; and furthermore, it can be used to define center of mass sections, in which the calculation of 
angular momentum gives an unambiguous intrinsic angular momentum. In particular the application of 
this intrinsic angular momentum to the two examples presented above, gives zero at any center of mass 
section. 

Incidentally, note also that for general spacetimes, inferring how much angular momentum is radiated 
by differentiating with respect to time eq. U]) will give results strongly affected by the issues discussed 
here. 
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